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Abstract. The key element in time-dependent density functional theory is the 
one-to-one correspondence between the one-particle density and the external 
potential. In most approaches this mapping is transformed into a certain type of 
Sturm-Liouville problem. Here we give conditions for existence and uniqueness 
of solutions and construct the weighted Sobolev space they lie in. As a result the 
class of ?;-representable densities is considerably widened with respect to previous 
work. 
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1. Introduction 

Imagine an interacting many-particle quantum system on a bounded space domain 
C K'^ within a time interval [0, T] governed by Schrodinger evolution from a fixed 
initial state ^'o G L^(O^), where TV is the number of particles. Using straightforward 
computational techniques one soon runs into complexity issues because of the multi- 
dimensional structure of the problem. One strategy to bypass this is given by time- 
dependent density functional theory (TDDFT). The key ingredient is to switch from 
the wave-function \1/ as a fundamental functional variable for evaluating obscrvablcs 
to the onc-particle density n : [0, T] x f2 K>o. That there exists indeed a one-to-one 
relation \I> o ri was first shown by Runge and Gross [1] for the time-dependent case 
but without considerations towards the domains of the involved mappings. 



In order to calculate the one-particle density n without the full complexity of 
the interacting Schrodinger problem one can rely on two possible schemes. Firstly, 
quantum fluid dynamical approximations like Thomas-Fermi theory and secondly, 
the usually adopted Kohn-Sham scheme. [2] Therein one substitutes the interacting 
quantum system by a non-interacting system in an external effective scalar potential 
i; : [0, T] X — 7> R leading to exactly the same one-particle density. This amounts to 
the fundamental question of (non-interacting) z;-representability. In order to guarantee 
the existence of such an effective potential one uses the divergence of the local force 
equation [J 

- V(nVv) = q[v] - dfn. (1) 
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Here the term q[v] : [0,T] x ft ^ R gives the divergence of the internal local 
forces, i.e. kinetic terms and all interactions. It is defined as an expectation value 
and thereby involves the state ^'[u] evolved under the influence of the 
external potential v. The dependence of q[v] on v is therefore highly non- linear. 

To solve equation ([T]) one will linearize its right hand side, e.g. by assuming 
analyticity in time of n and v hence finding a set of coupled linear partial differential 
equations. [3] A different approach currently under consideration by the authors tries 
to eliminate the analyticity restriction by pursuing a fixed point scheme to solve ([T]). 
[4] In any case one arrives at linear partial differential equations of Sturm-Liouville 
type for any time t £ [0, T] 

-V(nVv)=C. (2) 

The time-dependent generalization of density functional theory is nowadays a 
widely used technique in different fields of physics. [6] But in contrast to the ground- 
state theory there have not been a lot of rigorous mathematical investigations of 
TDDFT. As pointed out above the linear Sturm-Liouville equation is of fundamental 
importance to the foundations of the theory. Therefore it is crucial to have rigorous 
results concerning the properties of ^ in order to build a sound mathematical basis 
for TDDFT. So far uniqueness and existence of solutions w to ([2]) was shown in [5] 
but under the condition that the density is not only bounded but gapped away from 
zero for the whole domain il. In this work the authors considerably widen the class 
of allowed densities and give exact domains for the involved external potentials v as 
well as the inhomogeneity Further, we discuss implications on the involved external 
potentials and the one- particle densities within the different approaches [3j |4] to solve 
the non-linear equation ([T]). 

2. A problem-adapted weighted Sobolev space of potentials 

In the search for solutions to ([2|) we consider weak solutions, which are defined by 
adjoining an arbitrary u by means of the standard scalar product in all space- 
variables. 

-{u,VinVv)) ={u,0 (3) 

If we consider only potentials vanishing at the border of partial integration 
defines a bilinear form Q by 

Q{u, v) =^ (Vm, nWv) = {u, () ■ (4) 

The questions of existence and uniqueness of a solution w to can now be 
answered by the theorem of Lax-Milgram. [7] 

Theorem 1 Let Q be a coercive continuous bilinear form on a Hilbert space % . Then 
for every continuous linear functional C, on T-L, there exists a unique v £ % such that 
Q{u,v) = holds for all u£'H. 

A bilinear form Q is said to be coercive if there exists a constant c > such 
that Q{u,u) > c||u||^ for all u E Ti. Continuity means we find a C > such that 
Q{u,v) < C\\u\\-u ■ for all u,v & H. H as a Hilbert space should now be 

chosen in a way that makes Q coercive as well as continuous. The use of the scalar 
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product from before suggests H C Observing that (•, •) and Q can be naturally 

combined to a bilinear form 

(u, v)'H = (tt, v) + Q{u, v) = (u, v) + (Vm, nVv) (5) 

we can ask if this yields an adequate scalar product to define our Hilbert space H, i.e. 
we have to check if 

W^Wn = V {u, u)u = V + (^"^ "-^■") (6) 

is actually a norm. One obvious restriction is n > but this is certainly true for 
a one-particle density. The second restriction is such to make sense of the notation 
"Vm" if only in a distributional sense. Because of f2 bounded it holds L^{Vl) C L^(f2) 
and as the norm of % is constructed in such a way that ||w||2 < W^Wh we have the 
supposed property T-L C Remembering that elements of L\^^ can naturally be 

identified with distributions we are led to the following chain of inclusions allowing 
elements of Ti. to be distributionally differentiated. 

n c L'^in) c L\n) c Lioc(f^) c v'{n) (?) 

If n = 1 then || • is just the norm of the Sobolev space VF^'^(il) = H^{il) and 
written as || • ||i.2 so we adopt the notation H^(ri,n) for the complete normed space 
equipped with || • \\-u = \\ ■ \\i,2.n and call it a "weighted Sobolev space". [51 [H] 



Let us get our notation of the different norms involved straight by defining them 
all for general p £ [1, oo). Note the use of the weighting function n only in the second 
term of the definition of || • ||i,p,ri. 



\u\Pdx 



ri 


f \u\Pndx) 




\Jn J 



(8) 
(9) 
(10) 
(11) 



We still must not forget the restriction to functions which vanish at the border 
of il in order to justify the integration by parts used to derive This is of course 
not generally true for elements of H^(n,n) but can be met if one adopts the usual 
definition of Hq{V,) to weighted Sobolev spaces. In that we take the space of infinitely 
differentiable functions on fl with compact support C^(ri) and form the closure under 
our weighted Sobolev norm || • ||i,2,Tt. The resulting space H = (fi, n) equipped with 
scalar product ([5]) is complete and thus a fuU-fiedged Hilbert space of functions which 
vanish at the border of fJ. This will be the main space of our further investigations. 



Let us also define the dual of this Hilbert space H~^{n,n) — (iJp (il, n))', i.e. 
the space of linear continuous functionals on HQ(rt,n). In this we follow standard 
notation, cf. [5] 3.12 and 3.13. from Theorem [T] is thought of being an element of 
this space but we will rather concentrate on regular distributions as we try to find 
weak solutions to 
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3. Embedding theorems 

To prove coercivity of Q we follow a strategy largely outlined in [10]. The idea is to 
continuously embed the weighted Sobolev space into a non-weighted one and further 
into an ordinary space. 

Definition 1 Let V, W be Banach spaces with V C W. We say that V is continuously 
embedded in W and write V ^ W , if there is a constant c > such that for all v £ V 



\\v\\w < c\\v\\v 

We say that V is compactly embedded in W and write V 
bounded sequence in V has a subsequence converging in W . 



(12) 

W , if additionally every 



Lemma 1 Let be bounded, p > q > 1 and the weighting function n such that 



e L^{n) for s= then 



Proof. Using Holder's inequality with ^ 



p-i 
p 



1 wc derive 



< 



= \\\Vu\Pn\\l \\n- 



and thus 



(13) 

(14) 
^ (15) 

(16) 

Now we can easily establish the inclusion, considering that LP{il) C L'^{n) for 
bounded. □ 

The following Lemma is part of the Rellich-Kondrachov Theorem. [S] 
Lemma 2 Let be bounded then we have a compact embedding 

t¥™^«(ri) L'-in) 

for 1 > > - — ^ , provided m > 1 and mq < d. 



|Vu||g < c||Vu||j,,„. 



(17) 



Theorem 2 Let be bounded with dimension > 2 and the weighting function n such 
that 7i~^ S L^(r2) then we have a compact embedding 



(18) 



Proof. In the case of Lemma [T] choose p ~ 2,q ^ | which leads to s = 2. Now we can 
choose r = 2 and use Lemma [2] to get the following sequence of embeddings. 



(19) 



□ 
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4. Weak solutions to the Sturm-Liouville problem 

Theorem 3 Let Q be bounded with dimension > 2 and the weighting function n such 
that n^"^ G L^{Vl). The bilinear form Q defined on the Hilbert space HQ{Vl,n) by ^ 
is continuous as well as coercive and therefore admits the use of Theorem{J\ 

Proof. Part 1: Continuity. 

\Q{u,v)\ = |(Vw,nVi))| = KV^Vu, V^Vn)| (20) 

< ||^/^V^i||2 • \\^/T^Wv\\^ = \\^u\\^^^^ ■ \\^v\\^^^ (21) 

< Il"lli,2,.-Il«lk2,„ (22) 
Part 2: Coercivity. We need to show Q{u,u) > c||u||^2n f^i' ^ HQ{fl,n). 

We start by observing Q{u,u) = ||Vu||2 „■ First we use the result in the proof of 
Lemma [T] for the case p = 2, q = |. 

||V?/||4 <Ci||Vw|j2,n (23) 

As by Lemma[l]all elements of Hq (fl, n) are also in VFq ''(SI) we can apply the Poincare 
inequality (cf. |8] 6.30) \\u\\q < c||Vu||q to the l.h.s. It immediately follows 

< C2||Vli||2,„. (24) 

If we combine results ([23)) and ([24]) after taking (•)^ we get 

||^i||| + ||Vi.||! < (4 +c|) \\yu\\l^ (25) 

which leads us straight to the norm of W'^'^{n) and 

< C3||Vw||2,„. (26) 
Now we make use of Lemma [21 which tells us in the case m = l,q — |,r = 2 

ll"!|2 <C4||u||i,4 (27) 

and therefore the Hardy inequality |11] 

||w!|2 <C5||Vw||2,„ (28) 

holds. It is now easy to arrive at the desired inequality by squaring and adding 
another ||Vu||2_„. 

II«IIU«< (ci + l)!|V^i||l„ (29) 

□ 

Closely related we have the following Theorem as a solution to the general 
eigenvalue problem for Q. 

Theorem 4 Given a continuous and coercive bilinear form Q on H}f{Q,,n) under 
the conditions of Theorem [^ there is a monotone increasing sequence {\m)mefi of 
eigenvalues, 

< Ai < Aa < A,„ ™^ cx), (30) 

and an orthonormal basis {em\m&i C H^{il,n) of L'^lTt) such that for all u G 
Hl{Q,,n) and all m e N 

Q{u, era) = Am(w, e™). (31) 

Proof. By Theorem [2] we have a compact embedding Hl(yi,n) L^(fi) and this 

makes Theorem 6.3.4 in [7] applicable which yields just the given proposition. □ 
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5. Conclusions and Implications 

We conclude that it is possible to find a unique solution vt G HQ{^l,7^t) to ^ under 
the restrictions given in Theorem[3]and Q S H~^{p.,nt) for all t e [0,r]. Note that we 
made the time dependence of the involved sizes and thereby of the whole Hilbert space 
explicit by adding an index t. Even the space-domain f2 can be made time-dependent 
in this fashion. Further we have shown the existence of an eigenbasis of L^(ri) of the 
bilinear form Q from above with eigenvalues A,„ > 0. 

With respect to previous considerations [5] , where the linearization scheme of ([1]) 
poses restrictions only on the initial density no, we have extended the set of allowed 
densities. Especially we have lifted the restriction e < < M with M > e > for 
v-representable densities. However, in [4] one uses a linearization of ([T]) on the whole 
time interval. Hence we conclude that 

{{nt)tmT]\^t > o,nt e L\n),n^^ e L\n)} (32) 

is a sufficiently constricted set of densities to guarantee the existence of solutions to the 
linearized equation ^ for the whole time interval. This is an important requirement 
to make the fixed point approach presented in [3] rigorous and also affects previous 
considerations in this matter, e.g. [3]. 

Finally we want to discuss the condition ( e H~^{il, n). (We will drop the index 
t again in these last considerations.) Going back to ([T]) the inhomogeneity C of the 
linearized equations [3l |4] is noted to depend on q[v] ~ dfn. So in [U H] one will 
naturally assume that q[v] as well as dfn are in H^^{Q,n). The latter can be seen 
as an additional condition for w-representable densities. By virtue of the continuity 
equation dfn = —V • j [6] we transform this condition, which reads as {u, dfn) < oo 
for all u G Hq(Q,, n), into a more elementary form. 

{u,d^n) = -(w, V • dtj) = {Vu,dtj) - (^V^2\/u,n-htj) (33) 

The components of y/nVu are all elements of L^(51) because of u e HQ{n,n), so 
the question remains if this is also true for n^^dtj. This means we need to have a 
finite integral involving the force density dt'} 

lM!da;<c30. (34) 
n 

Note that this very much resembles the so called Weizsacker term from Thomas- 
Fermi theory which converges for finite kinetic energies. [12] 

y|Vy^|2dx = iy"^^d.T<oo. (35) 



Therefore (|34[) might prove useful in relating our condition on d^n to physical 
quantities as it is the case for the Weizsacker term. 
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